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The ring-diagram partial summation (or RPA) for the ground-state energy of the 
uniform electron gas (with the density parameter r s ) in its weak-correlation limit 
r s — > is revisited. It is studied, which treatment of the self-energy T,(k,uj) is 
in agreement with the Hugenholtz-van Hove (Luttinger-Ward) theorem fx — /zo = 
fj,) and which is not. The correlation part of the lhs has the RPA asymptotics 
alnr s + a' + 0{r s ) [in atomic units]. The use of renormalized RPA diagrams for 
the rhs yields the similar expression alnr s + a" + 0{r s ) with the sum rule a' = a" 
resulting from three sum rules for the components of a' and a". This includes in the 
second order of exchange the sum rule /X2 X = [P. Ziesche, Ann. Phys. (Leipzig), 
2006]. 
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I. INTRODUCTION 

Although not present in the Periodic Table the uniform or homogeneous electron gas (HEG) 
is still an important model system for electronic structure theory, cf. e.g. In its 

spin-unpolarized version, the HEG ground state is characterized by only one parameter r s , 
such that a sphere with the radius r s contains on average one electron j^jj]. It determines 
the Fermi wave number as k-p = l/(ar s ) in a.u. with a = [4/{9tt)} 1 / 3 » 0.521062 and 
it measures simultaneously both the interaction strength and the density such that high 



density corresponds 



on this limit cf. 



o weak interaction and hence weak correlation jjj. For recent papers 



Usually the total ground-state energy per particle is written as 



(here and in the following are wave numbers measured in units of kp and energies in kp) 

e = e + e x + e c , e = - • -, e x = -- • — , e c = {ar s ) 2 [a In r s + b + b 2x + 0(r s )}, (1.1) 

5 2 4 7r 

where e is the energy of the ideal Fermi g the exchange energy in lowest (1st) 

order, and e c is referred to as correlation energy given here in its weak-correlation limit 
with a = (1 - ln2)/7r 2 « 0.031091 after Macke [9] and b « -0.0711 after Gell-Mann and 
Brueckner n|. e c contains also the 2nd-order of exchange with b 2x ~ +0.02418 after 
Onsager, Mittag, and Stephen [llj. Notice that e = k F e = e/(ar s ) 2 gives the energy in a.u., 
e.g. the energy in zeroth order and the lowest-order exchange energy are eo = 3/(10 a 2 r 2 ) 
and e x = — 3/(47rar s ), respectively. 

Revisiting how Macke jjj, Gell-Mann/Brueckner Q], and Onsager/Mittag/Stephen Q, Q| 
derived e c in its weak-correlation limit, it is shown here, that and how an analogous procedure 
- also called RPA (= random phase approximation) - applies to the self-energy E(fc, u>). This 
latter quantity determines (i) the one-body Green's function G(k,u), from which follow 
the quasi-particle dispersion and damping and the momentum distribution n{k) [13]. It 
furthermore appears (ii) in the Galitskii-Migdal formula for the potential energy [14] (C + 
means the closing of the contour in the upper complex tu-plane), 

1 f „. q> f duj 
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/ d(k 3 ) f p: e lujS G(k,u)Z(k,u) , 5 > , (1.2) 



c+ 

which is related to the total energy e through the virial theorem [1 

d 



(iii) Besides, Ti(k,u>) appears in the Luttinger theorem Im E(l,/i) = in the 

Hugenholtz-van Hove theorem fi — fi — E(l,/z) ^?|, and in the Luttinger- Ward formula 
for the quasi-particle weight zf jlf] ]: 



1 



1 - E' 



E' = Re 



du 



;i.4) 



U) = fJ, 



So E(l,o;) is related to the chemical potential /x, which can be calculated from e according 
to the Seitz theorem 



/5 1 
At = I - - -r. s — | e 



d 

3 3 °dr S/ 

supposed e is known as a function of r s . Thus from Eq. (jl.lj) it follows for //: 



1 

2' 



/'x 



ar. s 



7T 



/'c 



ar. 



" Id - ( —-a + 6 + & 2x 



Similarly as in Eqs. and it is 

E(&, a;) = E x (fc) + E c (k, u), E x (/c) 



1 + 



k 2 



2k 



In 



k 



k 



0{r s ) 



ar 



;i.s) 



(1.6) 



— • (1-7) 
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Notice that E(fc, u;) in lowest order (of exchange) does not depend on to. In particular, it 
is S x (l) = — ar s /7r, thus /i x = S x (l). Similarly, in 2nd order of exchange the sum rule 
H2x = E 2x (l, \) holds [20]. With E 2x (l, \) = (ar s ) 2 C2 X it takes the form 6 2x = c 2x . The 
asymptotic behavior E c (l, jj) = (ar s ) 2 [a In r s + c + c 2x + 0(r s )} and the sum rule 



-a + b = c 
3 



1.8) 



are a must as a consequence of the Hugenholtz-van Hove theorem. But the question is: 
which partial summation of Feynman diagrams has to be used for E c and what has to 
be used for /i ? The obvious answer to the first question seems to be E c = E r + ■ • • (the 
subscript "r" means "ring diagram"). Symbolically written it is defined as E r = Go - {v r — vq) 
in terms of the Feynman-diagram building elements Gq, v , and Q(k,u) [ = polarization 
propagator in RPA], cf. Eqs. (|A.2[) . (|A.4|) . and Fig. [TJ v r is the effectively screened 
Coulomb repulsion following from v r — Vq + VqQv v , see Fig. [21 The Feynman dia- 
grams of E r are shown in Fig. H3 To what extent this (naive) ansatz has to be changed in 
a particular way (to answer also the second question) will be discussed at the end of Sec. III. 



Naively one should expect that in the weak-correlation limit the Coulomb repulsion e 2 /r 

i n 

3J can be treated as perturbation. But in the early theory of the HEG, Heisenberg [8J has 
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shown, that ordinary perturbation theory with e c = e 2 + e 3 + • • • and e n ~ (ar s ) n does not 
apply. Namely, in 2nd order, there is a direct term e2d and an exchange term e2 X , so that 
e2 = e2d + e2x- Whereas e2 X / (ar s ) 2 , cf. Fig. |U is a pure finite number fe 2x (not depending on 
r s ), the direct term e2d logarithmically diverges along the Fermi surface (i.e. for vanishing 
transition momenta q): e 2 d — ► hig for g — > 0. This failure of perturbation theory has 
been repaired by Macke |9| with an appropriate partial summation of higher-order terms 
e3 r , e^, ■ ■ ■ (the subscript "r" means "ring diagram") up to infinite order. This procedure 
replaces the logarithmic divergence for q — > by another logarithmic divergence, namely for 
r s — > 0, cf. Eq. (jl.lj) . This simultaneously "explains", why perturbation theory fails. The 
coefficient a, first found by Macke j^, has been confirmed later by Gell-Mann and Brueckner 
[lfl| . who in addition to the logarithmic term numerically calculated two contributions to 
the next term b, namely b r and 6 2 d- More precisely, instead of e r = e 2 d + e 3r + ■ ■ • (notice that 
e2r = e2d) they calculated a more easily doable approximation e° = + + • • • (which is 
sufficient in the weak-correlation limit) with the result e° = (ar s ) 2 [alnr s + b r + 0(r s )], so 
that e r = e° + Ae2d + O(rf) with Ae2d = e2d — e° d = {o.r s ) 2 b2A + 0{r^). In summary, 



e c = e r + e 2x + 0(rj 

= (e 2d + e 3r + ■ ■ •) + e 2x + 0{r 3 s ) 

= (e° d + 4 + ■ ■ •) + (e 2d - e° d ) + e 2x + 0(r 3 s ) 

= e° r + Ae 2d + e 2x + 0(r 3 s ) 

= (ar s ) 2 {[a\nr s + b r + 0(r s )} + b 2d + b 2x + 0(r s )} . (1.9) 

The result is Eq. (jl.lj) with b = b r + b 2 d- This procedure is revisited in Sec. 2 and then in 
Sec. 3 applied mutatis mutandis to the on-the-chemical-potential-shell self-energy E(l,/i), 
the rhs of the Hugenholtz-van Hove theorem. This is a contribution to the mathematics 
of the weakly-correlated (high-density) HEG). It concerns the HEG self-energy in RPA, 
extending and completing the paper p. 
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II. THE TOTAL ENERGY 



The Heisenberg-Macke story starts with the 2nd-order perturbation theory, e 2 = e 2 d + e 2x . 
Its components are the direct (d) term e 2 d (with qo^O) and the exchange (x) term e 2x : 

2-3 r d 3 q d 3 ^ d 3 k 2 P 



e 2 d = -(ar s y 



e 2x = +{ar s 



(27T) 



3><2o 



Q 1 



q-(k 1 + k 2 + q) 



d 3 q d 3 ^ d 3 k 2 



h,2 < 1, |fei,2 + g| > 1 , (2.1) 



h )2 < 1, |fci, 2 + g| > 1 • (2-2) 



(2tt) 5 J q 2 (fei + fc 2 + q) 2 q ■ (^ + fc 2 + q) 
P means the Cauchy principle value. (Notice the prefactor —1/2 and the replacement 
q A — > g 2 (fc! + fc 2 + q) 2 , when going from e 2( j to e 2x , and note that the 2nd-order vacuum 
diagram of Fig. El does not contribute.) As already mentioned, the integral ()2.2)1 has 
been ingeniously calculated by Onsager et al. [ill with the result e 2x = (ar s ) 2 6 2x , 6 2x = 
| In 2 — f^r ~ +0.0242. Unlike e 2x , the direct term e 2 d logarithmically diverges for q — > 0, 
i.e. along the Fermi surface. This is seen from 

2-3 



e 2 d 



d 3 q T( v 



(2-3) 



9><?0 



where the Pauli principle makes the function 

d 3 h d 3 k 2 p 



to linearly behave as I(q — > 0) 



q ■ (fei + fc 2 + q)' 

8tt 4 



fci2 < i, |fei,2 + q\ > i 



(2.4) 



a g + 0(g 3 ), see App. C, Eq. (RT2| . Thus e 2d 



(ar s ) 2 [alng 2 + •••], what agrees with (jl.lj) for q^ ~ ?%. The ring-diagram (or RPA) partial 
summation of Macke and Gell-Mann/Brueckner replaces the artificial (by hand) cut- 
off g by a natural cut-off q c ~ y / r s . This is made replacing the divergent direct term e 2 d by 
the non- divergent ring- diagram sum 

2 



16n 



' 2th ^2 n 



9c 



Q(q,v) 



'4ar, 



7T 



(2.5) 



For Q(q,r/), the polarization function in lowest order, is given in Eq. flA.4|) . With rj = iqu 
the contour integration along the real axis is turned to the imaginary axis: 



e r 



3 

8^ 



du J d(q 2 



q\ 



q 2 \nll + ^\R(q,u)-q 2 c R(q,u 



(2.6) 
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This has the advantage, that R(q,u) = Q(q,iqu) is a real function, being symmetric in u, 
cf. Eq. flA.2|) . Let us control Eq. ()2.6j) : The small-r s expansion of the u-integrand starts 
with (— l/2)(g c /g) 4 i? 2 (g, u), which just reproduces the 2nd-order direct term e 2 d with the 
help of the integral identity ()C.7j) . For r s — > 0, a direct numerical investigation of Eq. ()2.(i|) 
yields e r — > (ar s ) 2 (0. 031091 lnr s — 0.0711 + ■•■). This result is analytically rederived in the 
following. 



Namely, in the weak-correlation limit r s — > one can approximate R(q, u) ~ Q(qi— q)Ro(u) + 
■ ■ ■ with Ro(u) = 1 — w arctan 1/w, so that e r = e° + 0(r 3 ), where e° contains only the q- 
independent Rq(u) and its g-integration is restricted to < q < q\\ 



oo qi 



e r u = — / du \ dq 2 q 2 \n[q 2 + q 2 c R (u)] - q 2 \nq 2 - q 2 c R (u) 



3 

8^ 



o o 

oo 



J du -qtR%(u) 



In (%R Q (u) 



Oiri 



(2.7) 



(For a discussion of the divergent /convergent behavior of the g-series cf. jlOj . text after their 
Eq. (23).) With q 2 c = 4ar s /n it is 



3 f 

e° = (ar s ) 2 — / du R^{u) 
J 



4a 1 

In r s + In h In Ro (u) — 2 In q\ 

ix 2 



+ 0(r 



(2.8) 



So it results e° = (ar s ) 2 [a\nr s + b r — 2alngi + 0(r 3 )} with 



duR 2 (u) =a ^ 0.031091 



b r = a ( In 



4a 1 
~ ~ 2 

For the integrals cf. Eq. (jB.3|) . 



7T 



J du R 2 (u)\nR (u) « -0.045423 . (2.9) 



As it has been explained before and in Eq. ()1.9|) . the difference between the correct 2nd-order 
term of Eq. (|2.3|) and the first term in the expansion of e°, namely 

dq 
Q 



4d = -{ar s f 



2 - 3 J^-Jdu R 2 (u) + 0(r 3 s ) = -(ar s ) 2 2aj^- + 0(r 3 s 



7T' 



qo 



(2.10) 



gives 

Ae 2 d = e 2 d - e ! 



o 

2d 



-iar 



2 2-3 

7T a 



Leo 



dql{q) 
q 8nq 



1 <2l\ 



Vo l 



0(r, 2 



(2.11; 
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ar. 



I- 3 


j dq 


1 4vr 4 


J q 



/(g) 8n 4 



Q 



+ ] d j I M\ + 2a] d j\ + 0{r 



dq 
Q 



Ae 2 d = (aT s ) 2 [5 2d + 2a\nqi + 0(r s )} shows, that the sum e° + Ae 2 d does not depend on qi 
for r s — > 0. Besides, the first term of &2d is no longer divergent with go ~ > 0, therefore it can 
be set go = 0: 

i 



'2d 



dq 



4vr 4 J q 
o 

1 1 

4 + ^2 



/(g) 8tt 4 
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a 6(1 -g) 



11 8 . n 

— --ln2 + 2 In 2 2 

O 3 



-0.025677 . 



(2.12) 



Together it is b = b T + 5 2 d ~ —0.0711, what agrees with the above mentioned numerical 
evaluation of Eq. ()2.6|) . 



III. THE SELF-ENERGY 



Here - after the training of Sec. II - , it is aimed to calculate S c (l, fx) in the weak-correlation 
limit, where there is a scheme for £ c (l,/i) analog to Eq. (jl.9|) for e c with one difference. 
Namely, whereas the chemical-potential shift /i results from vacuum diagrams, the self- 
energy H(k,u) results from non-vacuum diagrams, which are functions of k and u, see the 
discussion at the end of this Section. 



In analogy to Eqs. (j2.1j) and ()2.2j) . the self-energy in 2nd order is £ 2 (&,a;) = Ti 2 d(k,uj) + 
£ 2x (/c,u;), the 2nd-order self-energy diagram of Fig. 03 vanishes. From (|A.6|) it follows for 
the direct term 

e(\e 1 + q\-l) 



J 2d 



1,W 



{ar s y 
2vr 4 



9><?0 



d 6 qd A k 2 



U-2 



+ 



- q ■ (ex + k 2 + q) + iS 
©(l-|ei + q|) 



(3.1) 
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Q(l - k 2 )Q(\k 2 + q\ - 1). 



^ - q ■ (d - k 2 ) - iS 

For the corresponding exchange term S 2x (l,u;) cf. Fig. H]and ref. [20J. where it has boon 
shown that £ 2x = Re S 2x (l, 5) = («r s ) 2 c 2x with the sum rule c 2x = b 2x ~ +0.0242. On the 
other hand, the direct term E 2 d diverges logarithmically for go — > 0. This is seen from 



J 2d 



Re S 2d ( 1 



1 



ar. 



2?H 



q>qo 



d 3 q ( , 
Y Ji ~ q) ' 



(3.2) 



8 



where the Pauli principle makes the function 

d 2 e x 



0(|ei + q[|-l)P e(l-|ei + g|)P 



q-(e 1 + k 2 + q) q- (e 1 - k 2 ) 



e(i- fe 2 )e(|fc 2 + 9 | - 1) , 

(3.3) 



to linearly behave as J(q — > 0) = 7c 3 a q + 0(q 3 ), see App. D, Eq. ()D.2|) . Thus S 2c j = 
(ar s ) 2 (a\nq 2 + • • •), what is for q 2 ~ r s in full agreement with the Hugenholtz-van Hove 
theorem (jl.4|) and the perturbation expansion of fi, which - because of (jl.5|) - gives fi 2 d — 
— {<yr s ) 2 (a\nr s + •••). In the ring-diagram partial summation the divergent direct term 
Y, 2 &(k,oj) is replaced by the non-divergent sum (its Feynman diagrams are shown in Fig. |3J) 

2 f d 3 q f drj Q(q, v) 



X 



q 2 J 2tti q 2 + q 2 Q(q,r]) 

«n _ |fc _u n \\ 1 

.(3.4) 



0(|fe + g|-l) + Q(l-\k + q\) 



_lo + 7] -\k 2 - q - (k + \q) +i5 uo + r] - \k 2 - q ■ (k + \q) - W_ 
Next, this expression is carefully controlled: 

(i) If the term q 2 Q(q, rj), which describes the RPA screening of the bare Coulomb interaction 
l/q 2 , is deleted, then E r (fc,u;) changes to £ 2 d(^,^), as it is seen from Eq. (|A.5|) . 

(ii) Use of Eq. ()3.4)1 in the Galitskii-Migdal formula (jl.2)) yields the ring-diagram summation 
for the potential energy, v T , which follows from e r through the virial theorem (jl.3j) . 

(hi) The expression ()3.4j) allows to calculate the derivative £[.(&, cu) = d'E r (k,uj)/duj. From 
S(. = Re £r(l> |) one obtains zp in RPA by means of the Luttinger-Ward formula (jl.4|) as 
Zp = 1 + + • • • with 



oc 



E' = / du ^4~r arctan ! + ... = -0.177038 r s + ■ ■ ■ . (3.5) 
7T 2 J Ro{u) u 







This is just the well-known RPA result for zp [24j. For the integral see Eq. (|B.iJ). 

After this controlling and training, E r = Re E r (l , ~) is derived from Eq. ()3.4|) in a similar 
way as e r in Eqs. (|2.7j) - (|2.12j) . The next steps again are substitution 77 = iqu and contour 
deformation from the real to the imaginary axis with x = e ■ e q and |e + q \ < 1 ± 5: 



7r 4 J q 2 J q 2 + q 2 R(q, u) {x + § ) - m 
(or.)" ?, u / ^ . 2(x + i) _ 



7T 4 J J q 2 q 2 + q 2 R(q,u) (x + § ) 2 + w 2 
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In the last step the u- and g-integrations are exchanged and it is used that R(q, u) is even 
in u, cf. Eq. flB.l|) ; so the imaginary part again vanishes. Next the (/-integration is specified 

as 



[ar s 



71' 



du I —- 



d 3 q R{q,u) 



oo oo 



[ar s ) 2 — J du J dq 



q 2 q 2 + q 2 R(q : u) 
R(q,u) 



dx 2{x + f ; 



2 (x + l) 2 + u 2 

1n (1 + 1) 2 + ^ 2 
(f -l) 2 + « 2 " 



(3.7) 



q 2 + q 2 R(q,u) 

Let us control Eq. (j3.7|) : The small-r s expansion of the ti-integrand starts with R(q,u)/q 2 , 
which just reproduces the 2nd-order direct term (j3.3J) with the help of the integral iden- 
tity (jD.7|) . In the limit r s — > 0, Eq. (j3.7|) numerically gives S r « (crr s ) 2 (0. 031091 In r s — 
0.081463 + ■ ■ ■). This result is analytically confirmed by the following. The asymptotic 
behavior for r s — > is determined by the lower integration limit q — > 0, therefore R(q,u) 
and In-- - can be approximated by Ro(u) = 1 — warctanl/w and L Q (u) = 2g/(l + u 2 ), 
respectively: 



err. 



7T" 



oo oo 

du J d{q 2 

n 



[R (u) + R^q 2 + ■ ■ -}[Lo(u) + L^q 2 + 



q 2 + q 2 [R (u) + R^q 2 



X? + 0(r s 3 ). 
(3.8) 



Finally the g-integration yields 



-(ar s ) — J du 
o 

oo 

2 [ 
-(ar s ) 2 — / du 

o 

2 °° 
-(ar s ) 2 — / du 

IT 



Ro(u) 



</l> 



dq 2 



1 + u 2 7 q 2 + q 2 R (u) 



R (u) 
1 + u 2 

Ro(u) 
1 + u 



Hq 2 + q 2 Mu)} I 

2 [2\nq 2 ~lnq 2 R (u)} 



(3.9) 



With Eq. (fB~3|) it turns out S° = (ar s ) 2 [a\nr s + c r - 2alng 2 + 0(r s )], 



, 4a 2 
a In 1 



7T 7T" 



^ o(M)ln t (M) - "0.035059 . 
1+u 2 



(3.10) 



The difference between the exact 2nd-order term of Eq. ()3.2|) and the first term in the 
g-expansion of namely 

2g 



^ = -{ar s ) 2 ± J^fdu Roiu) - 



qo o 



+ u 2 



-(or,)- — / — Trfl - In 2 



7T J « 

90 



(3.11) 
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yields 

AS 2d : 



J 2d 



J 2d 



or. 



2 


" oc 

/ 


dqJ(q) J 






q q 


go 


2 




jdq 


\J(q) 


7T 


3 


J q 

Lqo 


1 



dq 



tt(1 - ln2) 



(3.12) 



tt(1 - In 2) 



+ 



dq J (q) 
g g 



A£ 2d = (a^ s ) 2 [c2d + 2alng 2 + • • •] shows, that the sum S° + AE 2 d does not depend on q 2 
for r s — > 0. Besides the first term of c 2 d is no longer divergent with go — > 0, therefore it can 
be set go = 0: 



c 2 d = j / — 

7T 



c?g 
9 



£(9) 



- tt(1- In 2)9(1 -g) 



-0.046404 . 



(3.13) 



For the J(g)-integral cf. Eq. (|D.4jl . Together it is c = c r + c 2 d ~ —0.08146, to be compared 
with — |a + b = —0.08146. This is just the expected sum rule (|1.8|) . But the difference 
E r (l, g + A* x • • •) — S r (l, |) 7^ seems to disturb this sum rule c = — |a + b. This 'misfit' 
is removed by an additional partial summation replacing Fig. |3] by Fig. El i.e. replacing 
S r = G • (u r — v ) by = G x ■ (v T — Uq) with the renormalized one-body Green's function 

1 



\k 2 -E x (k)±i5 



10 



(3.14) 



see Fig. For k = 1 it is G x (l, w) = l/[w-|-// x ±i<y]. Thus E*(l, §+/j x +- • •) = S r (l, §)+• • • 
in the limit r s — > 0. So the conjectured relation (jl.8j) holds. This can be seen still in 
another way. 



Namely, note the similarities of Eqs. ()3.1()j) and (|2.9|) as well as Eqs. ()3.13|) and (|2.12jl . 
Their differences are 

1 2 
c r - b r = -a and c 2d - b 2 d = --a (3.15) 
o o 

using the identities ()B.5jl and (|D.8jl . respectively. Thus, with b = b T + b 2 & and c = c r + c 2 d, 
the sum rule ()1.8|) is proven once more. 



IV. SUMMARY 



Summarizing, the Hugenholtz-van Hove theorem /i — /xq = S(l, /i) takes for the HEG ground 
state in its weak- correlation limit r s — > the asymptotic form 

+ (ar,) 2 [alnr s + (-§a + b r + b 2d ) +b 2x + 0(r s )} = 
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- + (ar s ) 2 [a\nr s + (c r + c 2d ) +c 2x + 0(r,)] . (4.1) 

7T 

So the sum rules [with a = ^(1 — In 2) and 6 r , 6 2d , c r , c 2 d given in Eqs. (j2.9j) . (|2.12j) . (|3.1(JI) . 

1 2 

_a + 6 r = Cr , --a + &2d = c 2d , &2x = c 2x (4.2) 

hold. The last relation or /x 2x = S 2x has been shown in j^J. The sum rules ()4.2|) are 
relations between the Macke number a, the Gell-Mann/Brueckner numbers b T , fe 2d and the 
Onsager/Mittag/Stephen number 6 2x (which altogether desribe the r s — > asymptotics of 
the correlation energy e c ) on the one hand and corresponding numbers c r , c 2 d, c 2x of the 
on-shell self-energy E(l,/i) on the other hand. Eqs. ()4.1|) and ()4.2|) answer the question 
which partial summation of Feynman diagrams has to be used in the weak-correlation limit 
for the self-energy £ on the rhs of the Hugenholtz-van Hove theorem. They result from the 
renormalized ring-diagram (or RPA) partial summation (symbolically written as) £ £* = 
G x ■ v T with v T = vq/(1 — Qvq) and G x (k, uj) = renormalized one-body Green's function and 
Q(q,T)) = polarization propagator in RPA, see Eqs. (jH.14|) and (|A.4j) . respectively, and Figs, 
d d El CI (They not result from S w S HF = G ■ Vq, see Fig. |H1 as an alternative ansatz with 
G(k, lo) = full one-body Green's function of the interacting system (3].) Byproducts are the 
analytical representation of 6 2 d, a detailed description of the momentum-transfer or Macke 
function I(q) for the RPA vacuum diagrams (App. C), the introduction and discussion of an 
analog function J(q) for the RPA self-energy diagrams (App. D), and the proof of integral 
identities, which relate I(q) and J(q) to the polarization-propagator function R(q,u), cf. 
Apps. C and D. 
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APPENDIX A: ONE-BODY GREEN'S FUNCTION, PARTICLE-HOLE 
PROPAGATOR, AND 2ND-ORDER SELF-ENERGY 



In the following the identities (with z = x + iy) 



dz 




1 



2iri (z - zi){z - z 2 ) 



Z\ ~ z 2 

1 



for sign y\ = sign y 2 
for > and y 2 < 
for yi < and y 2 > 



(A.l) 



k Z 2 — Z\ 

for contour integrations in the complex z-plane are used (z = x + iy, C±= contour along the 
real axis, closed above or below with a half circle). The building elements of the Feynman 
diagrams are 



G (k,u) 



e(fc-l) 0(1 - k) 



cu-hk 2 + i5 u-hk 2 -i5 



5 > and vo(q) 



2" ' - - 2- - 1 

From Go follows the particle-hole propagator Q in RPA according to 

d 3 k r duj 



(A.2) 



(A.3) 



with the result 



Q{q,v) 



d 3 k 



-T7 + 



q(k + ~q) - 7] - i5 q(k + ~q) + 77 - 18 



e(l — fe)0(|Jfe + «| — 1). 



(A.4) 
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72|) and (|A.4|) used in the direct term of the 2nd-order off-shell self-energy 

2 r d 3 q f drj 



S 2d (A;,w) = (ar s ) 2 — y — J —Q(q,r})G (\k + q\,u + rj) 



g>go 



yields 



or. 



27T 4 7 Q 4 
<2><?0 



e(|fc + 9 |-i) 



w - |& 2 - q • (fe + k' + q) + i6 



(A.5) 



(A.6) 



+ 



e(i-\k + q \) 



uo-\k 2 -q.(k-k')- i6 



6(1 - k')Q(\k' + q\ - 1). 



This expression used in ()1.2|) yields t>2d = 2e 2 d in agreement with the virial theorem (J1.3|) . 



APPENDIX B: THE FUNCTION R(q,u) 



Q(q,f)) becomes real for imaginary rj: 



R(q,u) = Q(q,iqu) 



i + - 1 + ^ 



2q 

l+ q - 



(1 + 1)2+^2 
(f-l) 2 + M 2 



—u arctan 
V u 



1-3. 

— + arctan - 

u 



The function R(q, u) has the g-expansion R(q, u) = Rq(u) + q 2 Ri{u) + ■ ■ ■ with 

1 -5u 2 



Rq{u) = 1 — u arctan — , R\{u) 

u 

Here is a list of integrals: 



12(1 + u 2 ) 2 



R*(u) 



240(1 + u 2 ) A 



(B.l) 



(B-2) 



0?M Rq(u 







7T 

y 



du 



a w 0.321336, / lni2o(« 

o 

a w 0.482003, y rfw 



1 + M 



7T 
T 



Rn(u) In i?o(w 



oo 



i&M i 

— arctan — 
■u) it 



o 

oo 



+ U 2 



-3.353337 



-0.176945, 



-0.345751, 



R"(u) I 
du -£-4 arctan - « 4.581817 . 
Ro{u) u 



(B.3) 



(B.4) 



o v ' o 

The last but one integral appears in the weak- correlation limit of the quasi-particle weight 



2J. The identity 



7T° 



efax Rq(u) \iaR (u) 



1 3 

IT* - 2^ (u) 



(B.5) 



leads to the sum rule ()3.15|) . 
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APPENDIX C: THE FUNCTION I(q) 

Using cylindrical coordinates and the centre of the vector q as origin, Macke [l| succeeded 
to calculate I(q) explicitly as 

"29 fi ^ - 3 



I(q < 2) = tH 



ln2\ q - q - + i( 1 l + q -l + l) ]n (l + q - 

15 3 ) H 20 q\15 H 6 80/ V 2 

q I 15 y 6 80; V 2 



I(q > 2) 



7T 

30 



4 (22 + g 2 ) + -(q + 2) 3 (4 - 6g + q 2 ) In 1 + - 
Q \ U 

+ -(g-2) 3 (4 + 6g + g 2 )ln ( 1 - - 



5 



5, 



(C.l) 



Therefore /(g) is referred to as Macke function, cf. also |23J]. (The last two lines of (jC.lJl 
correct errors in |4j, Eq. (A.l).) /(g) has the properties 



I(q - 0) 

/(2) 
/"(2+) 



87T 2 7T 2 

— (l-ln2) g- — g 3 + ---, J(g->oo) 
3 6 



47T 2 

To" 



57T 



(13- 16 In 2) ^ 5.02598, J'(2) 
(2 + In 2) ss 14.1762, J"(2 - ) = 



8tt 2 



4tt\ 2 /1 2 1 
Ty I g^ + 5^ 



(-1 + 21n2) « -6.10012, 



2tt 2 



(7 -4 In 2) w -5.56305. (C.2) 



15 v 7 v y 15 

/(g) has a maximum of 7.12 at g = 1.36. /(g) and /'(g) are continuous at g = 2, but /"(g) 
has there a jump discontinuity of 27r 2 . This is because the topology changes from overlapping 
to non-overlapping Fermi spheres, when passing q = 2 from below. Its normalization is 



dq I (?) = -rr 
o 45 



-3 + 7r 2 + 61n2) « 19.3505. 



(C.3) 



/(g) is shown in Fig. 9. Multiplying the integral 



dq 



8tt 



7T 



— (l-ln2) gG(l-g)] = — [22 - 3?r 2 + 32 In 2 - 24(ln2) 2 ] 
3 9 



3.334856 
(C.4) 



with -3/(4tt 4 ) yields Eq. (l2~T2l . 



The original expression for /(g) arises from the diagram rules for e2d with Eq. (J3.3j) as 

/(g) = i(4vr) 2 Re /^Q 2 (g,??). (G.5) 
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One way is to insert (|A.4|) into (|C.5|) . It results (with Xi = q ■ {k i + \q),i = 1, 2) 

/ ( ,) = Re f ^!M% = / dyyp 

J x 1 + x 2 -iS J q ■ (fei + k 2 + q) 
in agreement with Eq. ()2.4|) . Another way is the analytical continuation and the deformation 
of the integration contour from the real to the imaginary axis with the advantage that 
R(q,u) = Q(q, iqu) is a real function. This yields 



I(q) = 2-Anq J du R 2 (q,u) 



(C.7) 



as an integral identity. 



APPENDIX D: THE FUNCTION J(q) 



Using again the method of Macke yields 



7T 



J(Q < 2) = jq 



|- 41n2 + |( 2 -|) + ln(l + f 



4n 

J(q>2) = —q[l 



1 ;i-g)(2 + g) 2 ln(l ' 



Q , 



J(q) has the properties 



7T 



J(q^O) = vr (l-l n2 )g-^g 3 + 
J(2) = y(l-ln2), J'(2+) = 



8tt 



l(l + g)(2-g) 2 ln ( 1 - - 
8g V q. 



. 4tt / 1 8 1 
J(g->oo =—-+—-+ 
3 \q z log 4 



7T 



-_(_! + 41n2), J'(T 



7T 



(D.l) 



-- -5 + 81n2 . 
o 



(D.2) 



(Notice I(q — > 0) = — J(g — > 0).) J(g) has a maximum of 1.3 at q = 1.9. J(q) is continuous 
3 

at q = 2, but t/'(g) has there a jump discontinuity. Its normalization is 



7T 



dqJ(q) = -(-3 + 7T 2 + 61n2). 
o 9 



(D.3) 



J(q) is shown in Fig. 10. Multiplying the integral 

°° , 

I -4 \J(q) -7r(l-ln2) qQ(l - q)] = -[10 - tt 2 + 8(1 - In 2) In 2] w 0.719405 (D.4) 
j o 8 
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with -2/vr 3 yields fjnSj) . 



The original expression for J(q) arises from the diagram rules for £2d with Eq. (jH.Hjl as 



J(q) = Re I d 2 e \ ^Q{Q,V) 



e(|e + g|-l)_ + G(l-\e + q\) 



q(e + \q) - r] - i8 q(e + \q) - r] + i5 



(D.5) 



One way is to insert (|A.4|) . It results [with x\ = q(ei + |q), x 2 = q{k 2 + \q)} 



J(q) 



d 2 e 1 
Air 



d 3 ko 



e(\ ei + q\-l)P , e(l-| ei + g|)P 



e(l-k 2 )e(\k 2 + q\-l) (D.6) 



X\ + X 2 X\ — x 2 

in agreement with Eq. (j3.4j) . Another way is the deformation of the integration contour 
from the real to the imaginary axis with rj = iqu. It yields 

u 2 + (l + §) 2 " 



J(q) = J du 
o 



In 



u 2 + (l-§) 2 



R(q,u) 



(D.7) 



as an integral identity. 



Comparing J(g) with /(g): 



57T 



7T" 



-a . 



(D.8) 



Note ^/(g -> 0) = J(q 0) = 7r 3 ag. The identity (ITXsll leads to the sum rule (ETToT) . 



Figures 




FIG. 1: Feynman diagrams for the one-body Green's function of the ideal Fermi gas Go(k,uj), the 
bare Coulomb repulsion vo(q), and the RPA polarization propagator Q(q,rj) as defined in Eqs. 
(|A~2l) -(TQl). 
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•lAMA/Wl? • • + - 




FIG. 2: Feynman diagrams for v r = vq + vqQv t , the screened Coulomb repulsion in RPA. 




FIG. 3: Feynman diagrams for the ring-diagram-summed self-energy S r = Go • (v T — Vq) as defined 
inEq. lEPl) . 




FIG. 5: Feynman diagrams, which do not contribute to e2 and S2, respectively. 





+ 



FIG. 6: Feynman diagrams for = G x • (v T — vq). For G x see Fig. |7|and Eq. (|3.14f) 
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FIG. 7: Feynman diagrams for the renormalized one-body Green's function G x = Go + GoT, x G : 
see Eq. (EHUl . 




FIG. 8: Feynman diagrams for E = Gvq with G = full one-body Green's function of the 
interacting system, G = G$ + GoT,G, S = full self-energy. The lowest-order term is S x = Got>o> 
therefore the correlation part is T,^ F = (G — Gq)vq. 




